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ABSTRACT 

The  numerical  solutions  are  obtained  for  skin  friction,  heat  transfer  to  the  wall  and 
growth  of  boundary  layer  along  the  flat  plate  by  employing  two  dimensional  Navier-Stokes 
equations  governing  the  hypersonic  flow  coupled  with  species  continuity  equations.  Flow  fields 
have  been  computed  along  the  flat  plate  in  CCF  atmosphere  in  the  presence  of  transpiration 
cooling  using  air  and  carbon  dioxide. 

1.  INTRODUCTION: 

The  main  focus  of  hypersonics  research  in  recent  times  is  on  the  re-entry  aerodynamics  in 
order  to  develop  protection  system  for  the  vehicles  from  severe  thermal  loads  encountered  during 
this  phase.  The  aerodynamic  heating  is  a  major  design  concern  for  vehicles  flying  at  hypersonic 
Mach  numbers.  The  thermal  load  problem  is  handled  by  using  an  effective  cooling  system.  The 
transpiration  cooling  in  which  a  cold  fluid  is  injected  into  the  boundary  layer  through  a  porous 
media  such  that  coolant  issues  forth  as  a  continuous  mass  is  an  important  cooling  technique. 
Many  experimental11 1  as  well  as  analytical  and  numerical  studies  of  transpiration  cooling  have 
been  reported  in  the  literature12,31.  Complete  Navier-Stokes  solutions  in  two-dimensional  form  for 
transpiration  cooling  in  earth's  atmosphere  using  air,  helium  and  carbon  dioxide  as  coolants  have 
been  reported  recently  using  CFD  techniques13'51.  The  focus  in  recent  times  has  been  on  the  planet 
Mars  and  many  probes  have  been  sent  to  explore  the  life  on  Mars.  This  paper  presents  an  analysis 
of  transpiration  cooling  in  carbon  dioxide  atmosphere  by  considering  a  simple  flat  plate  geometry 
flying  at  hypersonic  Mach  number  and  the  transpiration  cooling  is  achieved  by  injecting  air  and 
carbon  dioxide  gases.  The  N-S  equations  governing  the  hypersonic  flow  coupled  with  species 
continuity  equations  are  solved  numerically  using  finite  volume  technique.  The  numerical  results 
generated  using  this  CFD  code  have  been  verified  and  also  validated  by  conducting  experiments 
in  hypersonic  shock  tunnel  of  lISc. 

2.  GOVERNING  EQUATIONS  AND  NUMERICAL  ALGORITHM: 

Assuming  CO2  as  the  test  gas  flowing  at  hypersonic  Mach  number  while  the 
transpiration  cooling  is  achieved  by  injecting  a  coolant  gas  the  basic  Navier-Stokes  (N-S) 
equations  are  modified  with  appropriate  species  continuity  equations.  The  modified  non- 
dimensional  N-S  equations  in  the  conservative  integral  form  are  written  as: 

- 1 —  F  n  dS  =0  (1) 

dt  A  J 
s 

Where,  u  =  —  J  J  udA  is  the  average  value  of  the  column  vector  of  conserved  variables 

A 

u  =  [/).  pu,  pv,  pE,  /wi|  ]r  in  the  control  volume  and  F  is  the  flux  vector  defined  in  Ref.  [4],  Here,  p  is 
the  density,  u  and  v  are  the  velocity  components  in  the  x  and  y  directions,  E  is  the  total  energy  of 
the  mixture,  A  and  S  are  the  area  and  boundary  of  the  control  volume,  m/  is  the  mass  fraction  of 
the  species  land  n  is  the  unit  outward  normal.  The  viscosity  coefficient  p  for  the  gas  mixture  is 
expressed  as  a  function  of  temperature  based  on  either  Sutherland's  law  or  Wilkes  mixture  rule  [3]. 
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The  CFD  code  developed  is  based  on  the  finite  volume  method.  Discretising  and  applying  to  a 
cell  (i,  j)  the  above  equations  lead  to: 


dUj  |  1 

dt  Ay 


=  0 


(2) 


In  the  present  work,  Roe’s  scheme  is  used  to  determine  the  inviscid  fluxes.  The  quantities  such  as 
u,  v,  T,du/dx,  etc.,  are  required  on  an  interface  in  order  to  evaluate  the  viscous  fluxes.  In  the 


evaluation  of  the  derivatives  like  du/dxand  du/dy  at  an  interface  the  auxiliary  cell  approach  is  used. 


The  resulting  set  of  equations  is  treated  with  the  Lind-Jacobi  relaxation  algorithm. 

3.  RESULTS  AND  DISCUSSIONS: 

The  code  is  used  to  compute  flow  fields  at  Mx=6  &7  and  Re  =  106  and  5X  106.  The 
variation  of  normalized  skin  friction  coefficient  along  the  flat  plate  for  Mx=6  and  Re  =  106  for 
both  the  injecting  gases  for  different  blowing  parameter  is  shown  in  Fig  1.  The  C02  is  found  to  be 
more  effective  in  reducing  the  skin  friction.  The  variation  of  temperature  in  the  normal  direction 
and  the  variation  in  wall  temperature  for  the  same  Mach  number  and  Reynolds  number  are  shown 
in  Figs  2  &3.  The  results  clearly  show  that  the  thermal  boundary  layer  gets  affected  significantly 
due  to  blowing.  It  is  also  clear  that  the  injecting  gas  reduces  the  wall  temperature  significantly. 
The  experimental  investigations  have  also  been  carried  out  to  validate  the  CFD  results  and  will  be 
presented  in  the  full  length  paper. 
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Fig  1.  variation  of  skin  friction  for  injection  of  different  gases  for 
M  =  6.0  and  R  =  10  6 
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Fig2.  variation  of  temperature  profile  in  normal  direction  for 
M  =  6.0  andR  =10  6 
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Fig3.  Variation  in  wall  temperature  for  different  injecting  gases  for 
M  =  6.0  and  R  =  10  6 
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